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Abstract 

Perfect fluid Friedmann-Robertson- Walker quantum cosmological models for an 
arbitrary barotropic equation of state p = ap are constructed using Schutz's vari- 
ational formalism. In this approach the notion of time can be recovered. By su- 
perposition of stationary states, finite-norm wave-packet solutions to the Wheeler- 
DeWitt equation are found. The behaviour of the scale factor is studied by ap- 
plying the many-worlds and the ontological interpretations of quantum mechanics. 
Singularity-free models are obtained for a < 1. Accelerated expansion at present 
requires —1/3 > a > —1. 

PACS number(s): 04.20.Cv., 04.20.Me 

1 Introduction 

It is generally believed that the singularity which plagues the standard cosmological sce- 
nario may be avoided by taking into account quantum effects when, going backward in 
time, the Universe reaches the Planck scale. Since there is no consistent quantum theory 
of gravity until now, this possibility remains a speculation. However, important informa- 
tion concerning the singularity problem is expected to be obtained through the quantum 
cosmological approach In quantum cosmology the Hamiltonian formulation of gen- 
eral relativity is employed through the ADM decomposition of the geometry p|, and a 
Schrodinger-like equation, the Wheeler-DeWitt equation, is constructed which determines 
the wave function of the Universe as a whole. 

However, quantum cosmology suffers from many technical and conceptual problems. 
Technically, the Wheeler-DeWitt equation of quantum gravity is a functional differential 
equation defined in the so-called superspace, the space of all possible three dimensional 
spatial metrics, and no general solution in this superspace is known so far. In quantum 
cosmology this problem is circumvented by freezing out an infinite number of degrees 
of freedom by symmetry requirements, leaving only a few ones to be submitted to the 

^e-mail: flavio@cce.ufes.br 

^e-mail: fabris@cce.ufes.br 

•^e-mail: nivaldo@if.ufF.br 

"^e-mail: monerat@if.ufF.br 



1 



quantization process. This procedure defines the minisuperspace, where exact solutions 
can often be determined. On the other hand, upon applying the ADM decomposition 
general covariance is lost, and in most cases the notion of time in the quantum model 
disappears @]. Moreover, even if all these problems are solved, the question remains of 
the interpretation of the central object of quantum cosmology, the wave function of the 
Universe. Among the most popular interpretation schemes for the wave function of the 
Universe are the many-worlds ||^ and the ontological ^ interpretations of quantum 
mechanics. 

The many-worlds interpretation differs markedly from the Copenhaguen interpretation 
of quantum mechanics since, in a certain sense, the notion of probability is abandoned: 
All possibilities are actually realized and new universes are continously created by acts 
of observation according to the different possible eigenvalues obtained on measurement 
of an observable. But, in practice, the evolution of observables such as the scale factor 
is followed by means of the evaluation of expectation values. On the other hand, just 
like in the ordinary Copenhagen interpretation, the whole structure of Hilbert space and 
self-adjoint operators is kept intact. In its turn, the ontological interpretation of quantum 
mechanics makes use of the notion of trajectory. Thus, to a certain extent, the probabilis- 
tic concepts of the Copenhaguen school are also abandoned. This seems to be particularly 
interesting when we are discussing a quantum model for the Universe, since by definition 
the Universe encompasses everything, and the probabilistic concepts cannot be applied 
to it meaningfully. 

In a concrete study of a quantum cosmological model, a description for the matter 
content must be introduced. In principle, due to the quantum character of the problem, 
the matter content should be described by fundamental fields, as done in 0, for example. 
Predictions for the behaviour of the quantum Universe may be obtained, for example, by 
using the WKB approximation, establishing trajectories in phase space. However, general 
exact solutions are hard to find (even in the minisuperspace), the Hilbert space structure 
is obscure and it is a subtle matter to recover the notion of a semiclassical time |7|. 

In the present work, we describe matter as a perfect fiuid. This description is es- 
sentially semiclassical from the start, but it has the advantage of furnishing a variable, 
connected with the matter degrees of freedom, which can naturally be identified with 
time, leading to a well-defined Hilbert space structure. Another attractive feature of the 
phenomenological description is that it allows us to treat the barotropic equation of state 
p = ap with arbitrary a: General exact solutions can be obtained, constituting a nice 
"laboratory" for quantum cosmological models. 

To construct a quantum perfect fluid model, it is very convenient to use Schutz's for- 
malism 1^, ^ for a perfect fluid interacting with the gravitational fleld. In this formalism, 
certain velocity potentials are introduced giving dynamics to the fluid degrees of freedom. 
After some canonical transformations, at least one conjugate momentum associated to 
matter appears linearly in the action integral, and in this way a Schrodinger equation can 
be obtained with the matter variable playing the role of time. Therefore, all aparatus of 
ordinary quantum mechanics can, in principle, be employed in order to obtain predictions 
regarding the evolution of the Universe. 

Up to now, quantum perfect fluid models have been constructed only for the most 
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common equations of state, in particular those corresponding to dust, radiation and the 
vacuum [0, 0, |l^ . Predictions on the behaviour of the scale factor of the Universe 
have been made with the help of the many-worlds as well as the de Broglie-Bohm inter- 
pretations of quantum mechanics. For those special equations of state, universes with a 
bounce have been found, with absence of singularity; the classical behaviour is recovered 
for asymptotically large universes. 

Here, we generalize the previous investigations by studying quantum perfect fluid mod- 
els for any barotropic equation of state p = ap. The Wheeler-DeWitt equation is solved 
and wave packets are constructed. Next, using both the many- worlds and ontological 
frameworks, the behaviour of the scale factor is determined. For any value of a smaller 
than one a singularity-free bouncing Universe is obtained. Asymptotically, for large values 
of time, the classical behaviour is recovered. Although the results for the scale factor are 
independent of the interpretation scheme employed, the use of the ontological one allows 
us to verify explicitly that a repulsive quantum force appears as the Universe approaches 
the primordial singularity, leading to the bounce. The model predicts an accelerated ex- 
pansion today if —1/3 > a > —1. For a = 1 it is doubtful whether the quantum model 
exists at all, since we have been unable obtain finite-norm states due to divergences in 
the stationary wave functions. 

This paper is organized as follows. In the next section, the quantum cosmological 
model with a perfect fluid as the matter content is constructed in Schutz's formalism, and 
the Wheeler-DeWitt equation in minisuperspace is written down. The inner product and 
boundary conditions are given that insure self-adjointness of the Hamiltonian operator. 
In Section 3 the fiat case is considered. Wave packets are constructed and the expectation 
value for the scale factor is evaluated, as well as the Bohmian trajectories characteristic 
of the de Broglie-Bohm formalism. Sections 4 and 5 are dedicated to a brief discussion of 
the cases of positive or negative curvature. In Section 6 our conclusions are presented. 



2 The quantum model 

We need the Hamiltonian for a perfect fluid model in the formalism developed by Schutz. 
The starting point is the action for gravity plus perfect fluid, which in this formalism is 
written as 

S=[ + 2 f d^xVhhabK"^ + I d^x^^p , (1) 

JM JdM JM 

where K"''^ is the extrinsic curvature, and hat is the induced metric over the three- 
dimensional spatial hypersurface, which is the boundary dM of the four dimensional 
manifold M. Units are chosen such that the factor IGvrG becomes equal to one. The first 
two terms were first obtained in 0; the last term of (|l|) represents the matter contribu- 
tion to the total action, p being the pressure, which is linked to the energy density by the 
equation of state p = ap. In Schutz's formahsm p, ^ the fluid's four- velocity is expressed 
in terms of five potentials e, (, (3, 9 and S: 

U. = -{e,. + Cf3,u + eS,,) (2) 

/i 
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where /i is the specific enthalpy. The variable S is the specific entropy, while the potentials 
( and (3 are connected with rotation and are absent of models of the Friedmann-Robertson- 
Walker (FRW) type. The variables e and 6 have no clear physical meaning. The four- 
velocity is subject to the normalization condition 

f/"f/. = -1 . (3) 

The FRW metric 

ds'^ = -N'^dt^ + a?{t)aabdx''dx'' , (4) 

is now inserted in the action (j^). In this expression, N{t) is the lapse function and aah is 
the metric on the constant-curvature spatial section. Using the constraints for the fluid, 
and after some thermodynamical considerations, the flnal reduced action, where surface 



terms were discarded, takes the form |]TT 

r ' 2 O 

S= fdt -Q^ + 6kNa + N~'^''a^- ^(e + ^5)^/"+^ exp f - - 



(5) 



This reduced action may be further simplifled leading, by canonical methods to the 
super-Hamiltonian 

2 

7^ = -6fca + p"+^a-=^°e^ (6) 
24a 

where pa = —12aa/N and = —poU^Na^ , po being the rest mass density of the fluid. 
The canonical transformation 

T = -p^e- , Pt = P^V , e = e-(a + l)^ , p, = p, , (7) 

Pe 



which generalizes the one used in |Tl[, takes the super-Hamiltonian to the flnal form 

n = -^-6ka + ^ , (8) 
24a a3° ' ^ ^ 

where the momentum pt is the only remaining canonical variable associated with matter. 
It appears linearly in the super-Hamiltonian. The parameter k deflnes the curvature of 
the spatial section, taking the values 0, 1, —1 for a flat, positive-curvature or negative- 
curvature Universe, as usual. 

Imposing the standard quantization conditions on the canonical momenta and de- 
manding that the super-Hamiltonian operator annihilate the wave function, we are led to 
the following Wheeler-DeWitt equation in minisuperspace = 1): 

— - - lUka^^ + i24a^-=^"— = . (9) 

In this equation, t = — T corresponds to the time coordinate in a parametrization such 
that = a^", as follows from Hamilton's classical equations of motion |T^. By the way, 
the classical equation of motion for the scale factor is solved in a unifled form for any 
a G [0, 1] in terms of hypergeometric functions in [p^ ]. 
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Equation takes the form of a Schrodinger equation id'if /dt = H'^ . As discussed 
in Ijl^, |12|, in order for the Hamihonian operator H to be self-adjoint the inner product 
of any two wave functions $ and \1/ must take the form 



oo 

1— 3a/f,* 



$,^) = / a^-^°$*^c/a , (10) 







and restrictive boundary conditions must be imposed to the wave functions in the domain 
of the simplest ones being 



^(0,t) = or 



= . (11) 

a=0 



da 

The Wheeler-DeWitt equation can be solved by separation of variables. Indeed, 
writing 

*(a,t) = e-'^*e(a) (12) 

there results 

i" - UAka'^i + 24^0^"^"^ = , (13) 
where the prime means derivative with respect to a. 

3 The flat case {k = 0) 

As will be shown below, for a < 1 one can readily build wave-packet solutions to the 
time-dependent Wheeler-DeWitt equation (H) by superposition of stationary states. 

3.1 Stationary states 

For /c = the time-independent Wheeler-DeWitt equation (|13]) reduces to 

i" + 24Ea^-3°^ = . (14) 

If a < 1 it is possible to show that the parameter E is positive. Equation (0) admits a 
solution under the form of Bessel functions, leading to the following final expression for 
the stationary wave functions: 



-iEt r- T / v^96E 3(1-^) \ / 3{i-o) 

= e *^*v^ ci — -a 2 +C2F1 — ra 2 

3{i-<^) V3(l - a) / 3rr^V3fl-a) 



(15) 



The particular cases a = 0,1/3 and —1 have already been investigated in |TT|, 0, 



12, M. The above solutions are not valid for a = 1. In this special case equation (u_ 



becomes an Euler's type equation, the general solution of which takes the form 



V1-96-E V1-96-E 

cia 2 + 2 



(16) 



All of the above solutions must obey one of the boundary conditions (|TTp. The first 
one amounts to imposing C2 = 0, while the second one implies Ci = 0, except again for 
a = 1. The trouble with the two linearly independent stationary solutions for a = 1 is 
that, being just powers of a, their behaviour is irregular either at a = or at a = cxo. 
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3.2 Wave packets and the behaviour of the scale factor 

None of the stationary solutions found before has finite norm. Hence, wave packets must 
be constructed, by superposing those solutions, in order to obtain wave functions capable 
of describing physical states. The general structure of these superpositions is 

^(a,t)=/ A(E)'^Eia,t)dE . (17) 
Jo 

We specialize the discussion from now on to the case a < 1 and C2 = 0. Nothing of 
substance would be changed had we chosen ci = instead. Defining r = simple 
analytical expressions for the wavepacket are found if we choose the function A{E) to be 
a quasi-gaussian superposition factor: 

/"OO o q o/ \0 3(1 — a) 

^(a,t) = v^/ r^+V^'- *J,(ra^)rfr , (18) 
Jo 

where u = 7 is an arbitrary positive constant. The above integral is known 

||T6| , and the wave packet takes the form 

e 4B 

■^{a,t)=a 333^ , (19) 

where B = ■y — — aYt. Note that the norm of this wave function is finite only if 
a < 1. 

Now, we can verify what these quantum models predict for the behaviour of the scale 
factor of the Universe. In order to do this, we adopt first the many-worlds interpretation, 
and calculate the expectation value of the scale factor: 

^ /o°°ai-3"vl/(a,t)*avl/(a,t)rfa 



The above integrals are easily computed, leading to 

-9(1 -a)^ 



< a > (t) oc 



(32)272 



+ 1 



^ . (21) 



These solutions represent, for a < 1, a bouncing Universe, with no singularity, which goes 
asymptotically to the corresponding flat classical model, obtained from (|l|), when t ^ 00: 

a{t) oc t2/3(i-«) . (22) 
In order to fit observational evidence 0, 113 that the expansion of the universe is accel- 



erating, one must require — | > a > — 1. 

It is believed that the results obtained previously by means of the many-worlds in- 
terpretation scheme coincides with those that can be obtained using the ontological in- 
terpretation of quantum mechanics P, ^, in spite of a recent controversy on this issue 



p!9|, 20, 21 1 . In the ontological interpretation the wave function is written as 



6 



^ = Re 



iS 



(23) 



where R and S are real functions. Inserting this expression in the Wheeler-DeWitt equa- 
tion (pi), there resuhs 



dR 



+ 



1 



dS_ l_ 

'dt ^ 24ai-3" \daj 
dR dS 



— 1 +Q 



+ 



-R- 



dt 12ai-3« da da 24ai-3a Qa"^ 








(24) 
(25) 



where Q 
equation 



" 24ai-3" "^1^ quantum potential which corrects the Hamilton- Jacobi 

. When the quantum potential is more important than the classical potential, 



we can expect a behaviour deviating from the classical one. Notice that in the present 
case the classical potential is zero, since k = 0. 
The wave function (^) implies 



R 



S 



4y+(A)V-«)^t^ 

3 (1 - a)2a3(l-a)^ 



4-3a 
' 6Cl-a) 



a exp 



7a' 



3{l-a) 



4 



128 



(4 - 3a) 
+ :7Tz r arctan 



3 il-a)H 



3 1 



a) 



32 
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(26) 
(27) 



The Bohmian trajectories, which determine the behaviour of the scale factor, are given 
by 



Pa 



dS 
da 



Using the definition of pa-, taking the lapse function as = a' 
Bohmian trajectories becomes 



3o 



:28) 



the equation for the 



512- 



3(1 



t 



a] 



(29) 



which can be easily integrated to 

a{t) = ao 



7'+ 



{- 

V32 



'1 



(30) 



Oq being an integration constant. This is essentially the same behaviour found by compu- 
tation of the expectation value of the scale factor. The case a > 1 is classically forbidden 
since it predicts a speed of sound greater than the speed of light, and is also forbidden at 
the quantum level because the norm of the wave function becomes infinite and, strictly 
speaking, both the de Broglie-Bohm and many-worlds interpretations can not be applied. 



7 



It is curious, however, that formal Bohmian trajectories do exist for a > 1 and suggest a 
Universe that begins and ends in a singularity. 
The quantum potential takes the form 



Q(a,t) 



7 



1 — a 



32 



r 



(4 - 3a) 



37(1 — a)a 



3(l-a) 



(31) 



Inserting in this expression the trajectory 
terms of the scale factor only: 



the quantum potential can be written in 



Qia) 



7- 



32 



a 3(i_„)(4-3a) 



37(1 — a)aQ 



3(l-a) 



3(1-q) 



(32) 



From this expression, it is plain to see that for a < 1 the quantum effects become impor- 
tant near the bounce, while they become negligible for large values of a. Hence, asymp- 
totically the scale factor behaves classically. The force due to the quantum potential 
Fa = —dQ{a,t)/da is repulsive, leading to the avoidance of the singularity. 

The bad behaviour of the stiff-matter stationary solutions ([TBI) either at a = or at 
a = 00 has prevented us from finding finite-norm states by superposing them. This leads 
us to suspect that no perfect fluid quantum cosmological model exists for a = 1 and 
k = 0. 



4 The positive curvature case (k = 1) 

In this case, the quantum dynamics is governed by the Wheeler-DeWitt equation (||) 
with k = 1. We have been unable to find stationary solutions for arbitrary a, therefore 
we discuss separately the cases for which solutions could be found in terms of known 
functions. The case of radiation {a = 1/3) is omitted since it has already been treated in 

0- 



4.1 Cosmic strings {a = —1/3) 

Inspection of the Wheeler-DeWitt equation (|^) shows that although the geometry is closed 
the quantum dynamics is equivalent to that of the flat model. The stationary solutions 
take the form 



^(a,t) = e-*^*v^{Ci Ji/4(V-36 + QE a^) + ^2^1/4(^-36 + QE a^)} 



(33) 



A wave packet very similar to (|l^) with a = —1/3 can be constructed, and the behaviour 
of the scale factor follows the pattern of the fiat case. 
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4.2 Dust {a = 0) 



In this case the time-independent Wheeler-DeWitt equation (|T3D reduces to 



- C{a) + {-2AEa + UAa^) ^{a) = . 
In terms of the new variable x = 12a — E we find 



+ 



x' 

+ 



144 144 



e(a) = 0. 



(34) 



(35) 



Equation ( |35D is formally identical to the time- independent Schrodinger equation for a 
harmonic oscillator with unit mass and energy A: 



dx"^ 



+ 



-2A + w'^x^ 



^a) = 0, 



(36) 



where 2A = i?^/144 and w = 1/12. Inasmuch as the allowed values of A are n + 1/2, the 
possible values of E are 



n 



0,1,2,... 



En = Vl2(2n + 1) , 
Thus the stationary solutions are 

^„(a,t) =e-^^-Vn(12a-E„) 

where 



H„ 



X 



12/ 



\ -X.V24 



(37) 
(38) 

(39) 



with Hn the n-th Hermite polynomial. 

The wave functions (^Bp look like stationary quantum wormholes as defined by Hawk- 
ing and Page [Q. However, neither of the boundary conditions ( pH]) can be satisfied by the 
wormhole-like wave functions (|38D . Thus, at least in the dust case, our perfect fiuid model 
does not support static quantum wormholes, which are ruled out by the requirement that 
the Hamiltonian operator be self-adjoint. 



4.3 Stiff matter {a = 1) 

The general stationary solutions turn out to take the form 



(40) 



where K^, and I^, are modified Bessel fuynctions and u = — 96E/A. Since grows 
exponentially as a oo, here we must set C2 = and as a consequence the first of 
the boundary conditions (|TT]) is satisfied. Unfortunately, however, we have been unable 
to find explicit finite-norm solutions to the Wheeler-DeWitt equation by superposing 
stationary states because integrals over the order of modified Bessel functions are very 
hard to perform. 
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5 The negative curvature case (k = — 1) 



We managed to find stationary solutions to the Wheeler-DeWitt equation (P) with k = —1 
for the same values of a as in the positive-curvature case. 

5.1 Cosmic strings {a = —1/3) 

As in the positive-curvature model, in the present case although the geometry is open the 
quantum dynamics is equivalent to that of the flat model. A wave packet resembling ([T9|) 
with a = —1/3 can be readily constructed, the behaviour of the scale factor following the 
pattern of the flat case. 

5.2 Dust (a = 0) 

The stationary solutions are given by 



^{a, t) = e-^*(12a + Ey^^ \c,M^_^ , (^^^^^^) + C,W._^ , (^^^^^^) | (41) 

I 48 '4 i/ 48 '4 i/ J 

where M^^a and Wi^^\ are Whittaker functions, which are related to confluent hypergeo- 



metric functions [^]. The Whittaker functions in Eg . (PD do not automatically vanish 
at a = 0. Thus, in order to satisfy \l/(0,t) = it is necessary to take both Ci 7^ and 
C2 7^ 0, the same applying to the second of the boundary conditions (0). The difficulty 
in dealing with integrals over the order of Whittaker functions has prevented us from 
obtaining exphcit wave packets. 

5.3 Stiff matter {a = 1) 

The general stationary solutions turn out to take the form 

^(a,t) = e-'^'^[C^Jx{Qa^) + C^YxiQa")] (42) 



where A = a/1 + 96£'/4. If < -E < 1/96 any of the two boundary conditions (|TT|) can be 
implemented, but explicit wave packets could not be found by superposition of stationary 
states because very few results are known for integrals over the order of Bessel functions. 



6 Conclusions 

In this work we have investigated minisuperspace FRW quantum cosmological models 
with perfect fiuid for any value of the barotropic parameter a. In the case of fiat spacelike 
sections it has been shown that the models are completly solvable, except for the stiff 
matter case (a = 1), for which no quantum model could be constructed due to the 
divergent nature of the stationary wave functions either for large or small scale factor. 
The use of Schutz's formalism for perfect fluids allowed us to obtain a Schrodinger-like 
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Wheeler-DeWitt equation in which the only remaining matter degree of freedom plays 
the role of time. 

Superposing stationary wave functions, physically acceptable wave packets were con- 
structed in the case of flat spacelike sections. The time evolution of the scale factor has 
been determined in two different ways: Evaluating its expectation value, in the spirit of 
the many-worlds interpretation of quantum cosmology, and also computing the Bohmian 
trajectories of the ontological interpretation, in spite of controversies concerning the real- 



ity of such trajectories |2J]. In both cases the result is essentially the same for a < 1: A 
bouncing singularity-free Universe is obtained. The use of the ontological interpretation 
has allowed us to identify a quantum potential, from which a quantum force was com- 
puted. It acts repulsively as the Universe approaches the singularity, leading to avoidance 
of the singularity. In all cases, the classical behaviour has been recovered asymptotically. 
A Universe in accelerated expansion today requires —1/3 > a > —1. Near the bounce 
the behaviour of the scale factor is the same, irrespective of the value of a. 

The extension of the previous analysis to a non-flat spatial section leads to many 
technical new challenges. The case a = 1/3 has already been extensively studied in 



the literature |jTT], |T^, and we did not treat it here. In some other specific cases, the 
wave equation can be solved. For a = 0, the wave function is expressed in terms of 
Hermite polynomials {k = 1) or Whittaker functions {k = —1). Unfortunately, it seems 
unavoidable to perform numerical integration in order to obtain the wave packets from 
which the expectation value for the scale factor could be evaluated. The case a = — 1/3 (a 
cosmic string fluid) can also be solved, but it brings nothing new with respect to the flat 
case, since a cosmic string fluid mimics a curvature term. For a = 1, the same problems 
concerning the flniteness of the wave function appear, rendering the construction of the 
wave packet a very hard (or even impossible) task. 

We do not claim that our approach is superior to the more fundamental one based on 
quantum fields. In a subject still so far away from empirical tests, and so controversial 
as regards its physical meaning, we believe that several different approaches should be 
seriously considered and their consequences pushed as far as possible. This is what we 
tried to do here as concerns perfect fluid FRW quantum cosmological models. However, 
in order to construct more realistic models, as regards the earliest stages of the Universe, 
it is clearly important to consider fundamental flelds that can play crucial roles in that 
primordial phase. In particular, it may be interesting to include a conformal scalar fleld. 
We hope to present such a study in the future. 

Acknowledgements: We thank CNPq (Brazil) for partial flnancial support. 
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